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Abstract 

The general method for calculation of valence quark distributions in hadrons at 
intermediate x is presented. The imaginary part of virtual photon forward scattering 
amplitude on quark current with hadron quantum number is considered in the case, 
when initial and final virtualities of the current pf and ^-^e different, negative and 
large: | Pi |, | 1^ ^c^i where Rc is confinement radius. The operator product ex- 
pansion (OPE) in Pi,P2 up to dimension 6 operators is performed. Double dispersion 
representations in pf , of the amplitude in terms of physical states contributions are 
used. Equalling them to those calculated in QCD by OPE the desired sum rules for 
quark distributions in mesons are found. The double Borel transformations are ap- 
plied to the sum rules, killing non-diagonal transition terms, which deteriorated the 
accuracy in the previous calculations of quark distributions in nucleon. Leading order 
perturbative corrections are accounted. Valence quark distributions in pion, longitudi- 
nally and transversally polarized p-mesons and proton are calculated at intermediate x, 
0-2 ^ X < 0.7 and normalization points Q^ = 2 — 5 GeV^ with no fitting parameters. In 
cases of pion and proton the results are in agreement with found correspondingly from 
the data on the Drell-Yan process and deep inelastic scattering. Valence quark distri- 
butions in transversally and longitudinally polarized p-mesons are essentially different 
one from another and also differ from those in pion. 



1. Introduction 

Quark and gluon distributions in hadrons are not fully understood in QCD. QCD predicts 
the evolution of these distributions with in accord with the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) [l]-[3] equations, but not the initial values from which this evolution 
starts. The standard way of determination of quark and gluon distributions in nucleon is 
the following [4-8] (for the recent review see [9]). At some = Qq (usually, at low or 
intermediate ~ 2 — 5GeV^) the form of quark (valence and sea) and gluon distributions 
is assumed and characterized by the number of free parameters. Then, by using DGLAP 
equations, quark and gluon distributions are calculated at all and x and compared with 
the whole set of the data on deep inelastic lepton-nucleon scattering (sometimes also with 
prompt photon production, jets at high p± etc; for pion - with the data on Drell-Yan process). 
The best fit for the parameters is found and, therefore, quark and gluon distributions are 
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determined at all Q^, including their initial values q{Ql,x), g{Ql,x). Evidently, such an 
approach is not completely satisfactory from theoretical point of view - it would be desirable 
to determine the initial distribution directly from QCD. 

I present you based on QCD method of the calculation of u and d valence quark distri- 
butions in hadrons at low ~ 2 — 5 GeV^ and intermediate x. The idea of the method 
was suggested in [10] and developed in [11-13]. Recently, the method had been improved 
and valence quark distributions in pion [14] and transversally and longitudinally polarized 
p-meson [15] had been calculated, what was impossible in the initial version of the method. 
For valence quark distributions in proton the improved method [16] gave a much better 
agreement with experiment, than the earlier version [11]. The idea of the approach (in the 
improved version) is to consider the imaginary part (in s-channel) of a four-point correlator 
n(pi,P2, corresponding to the forward scattering of two quark currents, one of which 
has the quantum numbers of hadron of interest and the other is electromagnetic. It is sup- 
posed that virtualities of the photon g^, g'^ and hadron currents Pi,P2 are large and negative 
= Ig'^l ^ Ip^I, IP2I ^ Rc^7 where Rc is the confinement radius. It was shown in [11] that 
in this case the imaginary part in s-channel [s = (pi -|- q)^] of Il{pi,p2;qi,q') is dominated 
by a small distance contribution at intermediate x. (The standard notation is used: x is the 
Bjorken scaling variable, x — —q^/2v , v — piq). The proof of this statement, given in [11], is 
based on the fact that for the imaginary part of the forward scattering amplitude the position 
of the closest to zero singularity in momentum transfer is determined by the boundary of the 
Mandelstam spectral function and is given by the equation 

= -4-^/ (1) 
1 — X 

Therefore, if \p^\ is large and x is not small, then even at t = (the forward amplitude) the 
virtualities of intermediate states in t-channel are large enough for OPE to be applicable. 
So, in the mentioned above domain of q'^,g'^, and intermediate x I'mll{p,p2;q,q') 

can be calculated using the perturbation theory and the operator product expansion in both 
sets of variables q'^ — (p and pi, pi- The approach is inapplicable at small x and x close 
to 1. This can be understood for physical reasons. In deep inelastic scattering at large \q^\ 
the main interaction region in space-time is the light-cone domain and longitudinal distances 
along the light-cone are proportional to 1/a; and become large at small x [17,18]. For OPE 
vahdity it is necessary for these longitudinal distances along light-cone to be also small, that 
is not the case at small x. At 1 — x <S 1 another condition of applicability of the method is 
violated. The total energy square s = Q^(l/x— 1) + pi = —q"^ is not large at 1 — a; -C 1. 
Numerically, the typical values to be used below arc ~ 5GeV'^, pi ~ —IGeV"^. Then, 
even at x ~ 0.7, s ~ IGeV^, i.e., at such x we arc in the resonance region. So, one may 
expect beforehand, that our method could work only up to x ~ 0.7. The inapplicability of 
the method at small and large x manifests itself in the blow-up of higher order terms of OPE. 
More precise limits on the applicability domain in x will be found from the magnitude of 
these terms. 

The further procedure is common for QCD sum rules. On one hand the four-point cor- 
relator n(pi,p2; q, q') is calculated by perturbation theory and OPE. On the other hand, the 
double dispersion representation in pl,p\ in terms of physical states contributions is written 
for the same correlator and the contribution of the lowest state is extracted using the Borel 
transformaion. By equaling these two expression the desired quark distribution is found. 
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2. The method 



Consider the forward 4-point correlator: 

n(pi,P2;g,g') = -^ J d'xd'yd'ze'P^^+"-^ X {0\t{j\x), f\y), f\0), j\z)} |0) (2) 

Here pi and p2 are the initial and final momenta carried by hadronic current j^, q and 
q' = q + pi—p2 are the initial and final momenta carried by virtual photons. (Lorenz indeces 
are omitted). It will be very essential for us to consider non-equal pi,P2 and treat pl,pl as 
two independent variables. However, we may put = q'^ = q^ and t = [pi — P2Y = 0. We 
are interested in imaginary part of Il{pf,P2, q^, s) in s channel: 



ImU{p^,P2,q ^s) = — 



n(Pi,P2, Q ,s + i£)- H(pi,p2, q ,s-ie) 



(3) 



In order to construct representation of Imll{p'l,p2, g^, s) in terms of contributions of physical 
states, let us write for Imll{pl,pl, g^, s) the double dispersion relation in pl,pl: 



/mH(pl,p2,g^s) = a(g^s) + / '''^'f^ du+ [ 

J u — m J 



u — pI 



u-pl 



du 



+ / dU2 

J J {Ul — Pi)[U2 — P2I 



(4) 



The second and the third terms in the right-hand side (rhs) of (1) may be considered as 
subtraction terms to the last one - the properly double spectral representation. The first 
term in the rhs of is the subtraction term to the second and third ones. Therefore, (^) has 
the general form of the double spectral representation with one subtraction in both variables 
- Pi and P2- Apply the double Borel transformation in Pi,pl to @. This transformation kills 
three first terms in rhs of (HI) and we have 



00 00 

^Ml^Ml Imll{pl,pl, ^) = J / du2p{q^, s, Ui, U2)exp 



Ui U2 

Mf ~ Ml 



(5) 



The integration region over ui, U2 may be divided into 4 areas: 

I. ui < So; U2 < So 

II. Ml < So; U2 > So 

III. Ml > So; U2 < So 

IV. Mi,M2 > So 

Using the standard QCD sum rule model of hadronic spectrum and the hypothesis of 
quark-hadron duality, i.e. the model with one lowest resonance plus continuum, one may 
see, that area I corresponds to resonance contribution. Spectral density in this area can be 
written as 

p(mi, U2, X, Q"^) = gl ■ 27rF2(x, Q'^)5{ui - ml)6{u2 - ml), (6) 



where gh is defined as 



{0\jh\h)=gh 



(7) 
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(For simplicity we consider the case of the Lorenz scalar hadronic current.) If in 
Imll{pi,p2, q, q') the structure, proportional to P^Py [P^ = (^1+^2)^1/2] is considered, then in 
the lowest twist approximation F2{x, Q^) is the structure function, depending on the Bjorken 
scaling variable x and weakly on = —q^. 

In area (IV), where both variables -Ui 2 are far from resonance region, the non-perturbative 
effects may be neglected, and as usual in sum rules, the spectral function of hadron state is 
described by the bare loop spectral function p° in the same region 

p(mi,M2,x) = p°(ui,U2,a;) (8) 

In areas (11), (111) one of the variables is far from the resonance region, but other is in 
the resonance region, and the spectral function in this region is some unknown function 
p = il!{ui,U2,x), which corresponds to transitions like h — > continuum. The areas 11-111 
contributions are exponentially suppressed, and, using the standard hypothesis of quark- 
hadron duality, we may estimate them as a bare loop contribution in the same integration 
region. Equating physical and QCD representation of 11 and taking into account cancellation 
of appropriate parts in the left and right sides, one can write the following sum rules: 

-m^ ( 1 + 1 ) 

Im Il^QCD + Power correction = 27rF2(x, Q^)gle "1 

so so 

ImU'QCD = J j p\u^,U2,x)e-^-^ (9) 


It can be shown (see below), that for bare loop diagram V'('Wi; ^2, a;)) ~ 5(ui — U2), and, as a 
consequence, the areas 11 and 111 contributions are zero in our model of hadronic spectrum. 

It is worth mentioning that if we would consider the forward scattering amplitude from 
the beginning, put Pi = P2 = P and perform Borel transformation in p"^, then unlike (5), the 
contributions of the second and third terms in (^ would not vanish. They just correspond 
to the non-diagonal transition matrix elements and are proportional to 

{0\f\h*) {h*\f\x)f\0)\h) ^^{h\f\0) (10) 

p^ — mjf p'^ — 



From decomposition 



— m*j^ p^ — m\ m*j^ — m\ \p^ — m*^ p^ — m\ 



it is clear that in this case (PUj ) may contribute to the second (or third) term in (^) and 
after Borel transformation the contribution of the second term in ([TT|) has the same Borel 
exponent e~™'»''^^ as the lowest resonance contribution. The only difference is in pre-exponent 
factors: they are in front of the resonance term and Const, in front of non-diagonal 

terms. This difference was used in order to get rid of non-diagonal terms: application of 
the differential operator [d / d{l/ M'^)e'^h/^'^ to the sum rule kills the Borel non-suppressed 
nondiagonal terms, but deteriorates the accuracy and shrinks the applicability domain of the 
sum rule (particularly, the domain in x, where the sum rule is valid). 
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3. Quark distributions in pion 

It is enough to find the distribution of valence M-quark in tt"*", since d{x) = u{x). The 
most suitable hadronic current in this case is the axial current 

JmS = uif^lbd (12) 

In order to find the u-quark distribution, the electromagnetic current is chosen as u-quark 
current with the unit charge 

fj = ui,u (13) 
The bare loop Fig.l contribution is given by 

n-^- = - (2^ \ I ^ ik + pl-p,y ^^^'^ + - 

X Tr [^\k-ff,{k + q)^u{k+P2 - Piha{k ~ pi)^ (14) 
It is convenient to introduce 

P={pi+P2)/2, r = pi-p2,r^ = (15) 

The tensor structure, chosen to construct the sum rule is a structure proportional to 
P^PyPxPa/ V. The reasons are the following. The results of the QCD sum rules calculations 
are more reliable, if invariant amplitude at kinematical structure with maximal dimension is 
used. Different p\ ^ p2 are important only in denominators, where they allow one to separate 
the terms in dispersion relations. In numerators one may restrict oneself to consideration of 
terms proportional to 4- vector and ignore the terms ~ r^. Then the factor P^P^ provides 
the contribution of -^2(2;) structure function and the factor P\Pa corresponds to contribution 
of spin zero states. (The factor l/u is scaling factor : W2 = F2/v.) 
Let us use the notation 

n^vA. = {P^,p,PxP./vmplplx) + ... (16) 

Then ImT{{p\,p\,x) can be calculated from (|I4D and the result is [14]: 



3 00 00 ^, ^ . 

Imliip\,p\,x) = -x^(l - x) dui I du2 -, ^77 ^ (17) 

vr ' / i {ui-pj){u2-pl) 



The matrix element of the axial current between vacuum and pion state is well known 

{0\j,6\Tr)=tp,U (18) 



where /^r = 131Mey. The use of (^), (IT^), (P^ gives the sum rule for valence ti-quark 
distribution in pion in the bare loop approximation [14]: 

u.{x) = ^x(l - x)(l - e-^^"'y<l'''\ (19) 



where Sq is the continuum threshold. In ref. [14] the following corrections to (19) were ac- 
counted: 
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1. Leading order (LO) perturbative corrections, proportional to ln{Q'^/fi'^) , where /i^ is 
the normahzation point. In what follows the normalization point will be chosen to be equal 
to the Borel parameter /i^ = M^. 

2. Higher order terms of OPE. Among the latter, the dimension-4 correction, proportional 
to gluon condensate (0|^Gj]^ G^u\^) ^^^^ accounted, but it was found that its contri- 
bution to the sum rule vanishes after double borelization. There are two types of vacuum 
expectation values (v.e.v) of dimension 6: one, where only gluonic fields enter: 



TT 



and the other, proportional to four-quark operators 



(20) 



(21) 



It was shown in [14] that terms of the first type cancel in the sum rule and only terms of the 
second type survive. For the latter one may use the factorization hypothesis which reduces 
all the terms of this type to the square of quark condensate. 

A remark is in order here. As was mentioneed in the Introduction, the present approach 
is invalid at small and large x. No-loop 4-quark condensate contributions are proportional to 
(5(1 — x) and cannot be accounted. In the same way, the diagrams, which can be considered 
as a radiative corrections to those, proportional to S{1 — x), must be also omitted. 

All dimension-6 power corrections to the sum rule were calculated in ref.'s 14,15 and the 
final result is given by: 



3 M^^</M^ 



XUTriX] 



27r2 



J TT 



xHi 



X 



1 + 4a;/?T.(l 



X) 



X 1 — X 




a,{M^)-ln{Ql/M' 



2(1 - 2x)lnx 



(27r)4 . 37 ■ 26 ■ M6 x^{l-x)\ 
where uj{x) is 4-order polynomial in x, (for its explicit form see [14]), 

a = -(27r)2(0|V^^/'|0) 



(22) 



(23) 



Ut,{x) may be used as an initial condition at = Ql for solution of QCD evolution equations 
(DGLAP equations). 

In numerical calculations we choose: the effective LO QCD parameter Ag^^ = 
250MeV, Ql = 2GeV^, a^a'^ilGeV^) = OMGeV^. The value of the_latter is taken from the 
recent best fit of hadronic r-decay data [19] and corresponds to agiipip)"^ = (2.2 ± 0.7) ■ 10"^ 
GeV^. The continuum threshold was varied in the interval 0.8 < Sq < 1.2GeV^ and it was 
found, that the results only slightly depend on it. The analysis of the sum rule ( |2^ ) shows, 
that it is fulfilled in the region 0.15 < x < 0.7; the power corrections are less than 30% and 
the continuum contribution is small (< 25%). The stability in the Borel mass parameter 
dependence in the region OAGeV"^ < < O.GGeV^ is good. The result of our calculation 
of valence distribution in pion xUt^Ix, Q^) is shown in Fig.2. 



6 



Suppose, that at small x < 0.15 Ut,{x) ~ according to Regge behaviour and at 

large x > 0.7 ~ (1 — x)^ according to quark counting rules. Then, matching these 

functions with (^2]) , one may find the numerical values of the first and second moments of 
M-quark distribution 

1 

A^i = y u^{x)dx 1.1 (24) 



1 

M2 = J xu^{x)dx ^ 0.24 (25) 



Ail has the meaning of the number of w-quarks in 7r+ and should he Ai = 1. The deviation of 
(P^) from 1 characterizes the accuracy of our calculation. A^2 has the meaning of the part of 
pion momentum carried by valence w-quark. Therefore, valence u and d quarks are carrying 
about 40% of the total momentum. In Fig.2 are plotted also the valence w-quark distribution 
found in [4b] by fitting the data on production of fi~^fj^~ and e~^e~ pairs in pion-nucleon 
collisions (Drell-Yan process). 

4. Quark distributions in /o-meson 

The choice of hadronic current is evident 



f, = u7,d (26) 
The matrix element {p"^ | | 0) is given by 

2 

771 

{P^ I I 0) = ^e, (27) 
9p 

where mp is the p-meson mass, gp is the p — 7 coupling constant, gp/4:n = 1.27, is the 
p-meson polarization vector. Consider separately w-quark distributions in longitudinally and 
transversally polarized p-mesons. It was shown in [15], that for determination of w-quark 
distribution in longitudinal p-meson the most suitable tensor structure is that, proportional 
to P^PyPaPx, while M-quark distribution in transverse p can be found by considering the 
invariant function at the structure — P^Py^Ao-. In case of longitudinal p-meson the tensor 
structure, which is separated is the same, as in the case of pion. Since at = bare loop 
contributions for vector and axial hadronic currents are equal, the only difference from the 
pion case is in the normalization. It can be shown, that w-quark distribution in longitudinal 
p-meson can be found from (22) by substitutions — > mp, — > 7np/gp and therefore 

xu^p{x) = Ilgy-^h<y^^'xu^{x), (28) 

where xuTr{x) is given by (22). (The numerical values of and sq differ in the cases of 
pion and longitudinal p-meson). After separation of the invariant function at the structure 
P^iPu^Xa in box diagram we find w-quark distribution in transversally polarized p-meson in 
bare loop approximation: 



x(l — x) 

2 ^ ' 



^Mx) (29) 
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u 



Tl 



x)o has a minimum at a; = |, unlike Ut^{x)q and Up{x)Q. Another essential difference in 
comparison with pion and longitudinal p arises due to the fact that gluon condensate and 
(G^) operator contributions are nonvanishing in the sum rule for u^{x). 

Valence w-quark distribution in transversally polarized p-meson is given by 



xu. 



87r2 



ml 



X I (fo{x)Ei 





+4Zn(l - x) 



2{1 - 2x + 4:x'^)lnx^ 
ipo{x) 



^ (0|(a,/^)G2|0) 
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'(4a;- 1) 



M^x^ 



- 2« ■ 3^M^x3(l - ,)3 (0|^^/°"^::^^^A^A.|0)e(^) 



X) 



(30) 



where ^{x) and x(a;) are polynomials in x [15]. In numerical calculations the values of 
parameters were chosen: = 1 GeV^ (this value is in the middle of stability interval in 
M2), So = 1.5 GeV^, ^ 4 GeV^, (0 | (a./Tr)^^ | 0) = 0.006 GeV^ [20], (G^) vacuum 
expectation value was estimated according to instanton model [21]. The effective instanton 
radius, was taken to be equal Pc = 0.5 fm. The other parameters arc the same as in calculation 
of xUn{x). On Fig.3 are plotted Up{x),u^{x) and, for comparison, Ut,{x). As is seen from 
the Figure, the shape of u^ix) is quite different from u.,r{x) and Up{x) and indicates, that 
perhaps, u^{x) has two humps. Up{x) also has an interesting feature: the second moment 
of M-quark distribution is equal to 0.4. It means, that the part of p^ momenta, carried by 
valence - u and d quarks comprises 80% and for gluons and sea quarks is left about 20% only 
- a very unusual situation ! 

The main physical conclusion is: the quark distributions in pion and p-meson have not to 
much in common. The specific properties of pion, as a Goldstone boson manifest themselves 
in different quark distributions in comparison with p. 

5. Quark distributions in proton 

Consider the 4-current correlator which corresponds to the virtual photon scattering on 
the quark current with quantum number of proton: 

T'"'{puP2: q, q') ^-ij d'xd'yd'z ■ e^(fi-+"^) . {Q\T{r^{x), r/{y),rAQ), f){z)m, (31) 



where rj{x) is the three-quark current rj{x) = e {u°-C^\ii )'^^'-)\d^ [22], = 'U7y^M, = d'-^^d. 
As was shown in [11] and generalized to nonequal pi,P2 in [16] the sum rules should be written 
for invariant amplitude at the structure pp^Pv (it what follows it will be denoted by ImT/v). 
The bare loop diagram is shown on Fig.4. The result of its calculation after the double Borel 
transformation are the same as in the case for equal pi = p2 [11]: 



i(d) 



327r3 



where 



ip^(x) = x{l - x)^{l + 8x), ip^ix) = x{l - xf{l - 2x), 



(32) 
(33) 
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So is the continuum threshold 



E2{z) = 1 - (1 + ^ + z''/2)e~' (34) 
The sum rules for u and d-quark distributions q{x)'"'^'^^ in bare loop approximation are: 

xqix)^' = W W ■ ^2(-^) (35) 

— 2 

Here Xn is the coupling constant of proton with the current (0|?7|]9, r) = Xnv'''{p)^ = 
327r^A^. The first and second moments of quark distributions, with account of proton mass 
sum rule in the same approximation, coincide with quark model results, as it should be. 
QCD evolution was accounted as leading order (LO) correction to bare loop formulae (see 
[16]). In OPE the contributions of = 4 operator (gluon condensate) and d = Q operator 
(cts times quark condensate square) was calculated. I present here results as the ratio to bare 
loop contributions: 

u{x)i^G^) ^ {{aJ'n)G'^) vr^ (11 + 4x - 31a:^) _ so/M^\/p (_^\ (oa\ 
Uo{x) 12 a;(l-a;)2(l + 8a;) ^ 'W^^ ^ ^ 

(l_e-WM^)/i^,(^) (37) 

(38) 



do{x) M4 6 x^{l-x)^{l + 2xy 

^W)<xsmf _ 0LsC?{2\^ - 867a; + 172a;2 + 288(1 - x)ln2) 1 
Uq{x) ~ M6-8l7r-8x(l-x)3(l + 8x) E^is^jW^) 

d{x)as{i^^l^y _ Oistt^ilQ - 43a; + 36x'^) 1 



(39) 



do(x) M68l7ra;(l - x)3(l + 8a;) E2(so/M'^) 

On Fig. 5 are plotted the obtained u and (i-quark distributions at Qq = 5 GeV^ in comparison 
with those, found in [4a] from the data on deep inelastic lepton-nucleon scattering. The same 
values of QCD parameters (condensates) were taken as in cases of quark distributions in tt 
and p. The Borel parameter was chosen in the middle of stability interval, = 1.1 
GeV^. 

The theoretical analysis of the obtained valence quark distributions at = 5GeV^ 
showed, that -u-quark distribution is reliable at 0.15 < a; < 0.65, its accuracy is about 25% in 
the middle of this interval and decreases to 50% at the ends of interval; rf-quark distribution 
is reliable at 0.25 < x < 0.55 with an accuracy of about 30% in the middle and is given by a 
factor of 2 at the ends of interval. In the limit of these accuracies the theoretically calculated 
valence quark distributions are in agreement with those found from deep inelastic scattering 
and other hard processes data. 

The research described in this publication was made possible in part by Award No.RP2- 
2247 of the US Civilian and Development Foundation for the Independent States of Former 
Soviet Union (CRDF), by INTAS Grant 2000, Project 587 and by the Russian Found of Basic 
Research, Grant No. 00-02-17808. 
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Figure 1: Diagrams, corresponding to the 
unit operator contribution. Dashed hues with 
arrows correspond to the photon, thick sohd 
- to pion or rho current. 



Figure 4: Bare loop diagram in case of proton. 



xu(x) 



0.1 0.2 0.3 



xd(x) xu(x) 

Figure 2: Quark distribution function in pion 
(dashed hne) and the fit from [4b] - solid fine. 
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Figure 3: Quark distributions in longitudi- 
nal p, (sohd hne with squares), transversal p 
(solid line) and pion (dashed line) 



X X 
a b 

Figure 5: Valence u- and rf-quark distribu- 
tions in proton (dashed lines) in comparison 
with the fit [4a] of deep inelastic scattering 
data (solid hues) 
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